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Abstract: We consider the problem of reconstructing binary images from their horizontal and 
vertical projections. It is known that the projections do not necessarily determine the image 
uniquely. In a previous paper it was shown that the symmetric difference between two solu- 
tions (binary images that satisfy the projections) is at most 4a\/2iV. Here is the sum of the 
projections in one direction (i.e. the size of the image) and a is a parameter depending on the 

^ projections. In this paper we give a lower bound: for each set of projections that has at least two 

solutions, we construct two solutions that have a symmetric difference of at least 2a + 2. We also 
T— I show that this is the best possible. 

> 

o 

^ 1 Introduction 
o 

On 

O An important problem in discrete tomography is to reconstruct a binary image on a lattice 

j> from given projections in lattice directions [HI E] . Each point of a binary image has a value 

kj> equal to zero or one. The line sum of a line through the image is the sum of the values of 

the points on this line. The projection of the image in a certain direction consists of all the 
^ line sums of the lines through the image in this direction. 

For any set of more than two directions, the problem of reconstructing a binary image 
from its projections in those directions is NP-complete [S]. For exactly two directions, the 
horizontal and vertical ones, say, it is possible to reconstruct an image in polynomial time. 
Already in 1957, Ryser described an algorithm to do so [8]. He also characterised the set of 
projections that correspond to a unique binary image. Suppose F is uniquely determined 
and has row sums ri, r2, . . . , r^- For each j with 1 < j < maxj rj we can count the number 
: > j} of row sums that are at least j. Then these numbers are exactly the non-zero 
column sums of F (in some order). See also [HI Theorem 1.7]). 

Alpers et al. [H [2] studied the possible difference between a uniquely determined image 
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and a second image with almost the same projections as the first one. Their results were 
generalised in [3] and the same ideas were used to study the difference between two solutions 
of the same set of projections in We give an overview of the main theorems here. 

Consider given row sums TZ = (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„), and 
assume that there exists at least one binary image with exactly these line sums. Define V = 
(f 1, V2, . . . , Vn) as Vj = : r/ > j} for I < j < n. Let Fi the uniquely determined binary 
image with row sums TZ and column sums V. Let = ^™ ^ = X]j=i^i- Furthermore 
define the integer 

1 " 

a{n,c) = 

The parameter a indicates how close the line sums {TZ, C) are to line sums that uniquely 
determine an image. In particular, a = if and only if there is exactly one binary image 
with line sums {TZ,C). Intuitively, the larger a, the more possibilities there are for images 
that satisfy the line sums. 

Alpers et al. [U |2] proved that if F2 is an image with line sums (Jl,C) and a(7l,C) = 1, 
then the size of the symmetric difference between F2 and the uniquely determined set Fi 
is bounded: 

|Fi AF2I < V8N + 1-1. 

In [3] we generalised this to larger values of a(7l, C). For an image F2 with line sums (JZ, C), 
write a = a{TZ,C) and let p = |Fi fl F2I. Then 

|Fi A F2I < 2a + 2(a + p) \og{a + p) 

and 

\Fi A F2I < aVSN + 1 - a. 

The first bound is asymptotically sharp when a is large compared to p. The second bound 
is better when a is small compared to p. We have a family of examples in which this bound 
is achieved up to a factor pa- 
using these results, we can also consider the difference between two images with the same 
line sums. This must be bounded by twice the upper bound for the difference between one 
of these images and a uniquely determined image. Hence we get the following result from 

HI- 

Suppose F2 and F3 are two images with the same line sums (7Z, C). Write a = {7Z, C). Then 

IF2 A F3I < 2aV8N + 1 - 2a. 
As with the previous bound, there are examples for which this bound is only off by a factor 
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In this paper we consider the complementary problem: find the best lower bound for the 
symmetric difference between two solutions that you can at least achieve given a set of 
projections? For each set of projections that has at least two solutions, we construct two 
solutions that have a symmetric difference of at least 2a + 2. We also show that this bound 
is sharp. 



2 Definitions and notation 

Let -F be a finite subset of 1? with characteristic function (That is, ^(s;; y) = 1 if (2;, G 
F and x{x, y) = otherwise.) For i G Z, we define row i as the set {{x,y) E Z"^ : x = i}. We 
call i the index of the row. For j G Z, we define column j as the set {{x,y) G : y = j}. 
We call j the index of the column. Following matrix notation, we use row numbers that 
increase when going downwards and column numbers that increase when going to the right. 

The row sum ri is the number of elements of F in row that is = YljezXih j)- The 
column sum Cj of F is the number of elements of F in column j, that is Cj = Yliez xihj)- 
We refer to both row and column sums as the line sums of F. We will usually only 
consider finite sequences 71 = (ri, r2, . . . , r^) and C = (ci, C2, . . . , c„) of row and column 
sums that contain all the nonzero line sums. We may assume without loss of generality 
that ri > r2 > . . . > and Ci > C2 > . . . > c„. 

Given sequences of integers 71 = (ri, r2, . . . , r^) and C = (ci, C2, . . . , c„), we say that {TZ, C) 
is consistent if there exists a set F with row sums 7Z and column sums C. We say that the 
line sums (7Z, C) uniquely determine such a set F if the following property holds: if F' is 
another subset of 7? with line sums {7Z,C), then F' = F. In this case we call F uniquely 
determined. 

We will now define a uniquely determined neighbour of a set F. This is a uniquely deter- 
mined set that is in some sense the closest to F. See also [31 Section 4]. 

Definition 1. Suppose F has row sums > r2 > . . . > and column sums Ci > C2 > 
. . . > c„. For I < j < n, let vj = #{/ : ri > j}. Then the row sums ri, r2, . . . , Vm and 
column sums Vi, V2, ■ ■ ■ , f„ uniquely determine a set Fi, which we will call the uniquely 
determined neighbour of F. 

Note that if F' is another set with row sums ri, r2, . . . , and column sums Ci, C2, . . . , Cn, 
then Fi is a uniquely determined neighbour of F' if and only if it is a uniquely determined 
neighbour of F. Hence Fi only depends on the row and column sums and not on the 
choice of the set F. We will therefore also speak about the uniquely determined neighbour 
corresponding to the line sums {7Z,C), without mentioning the set F. 
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Suppose line sums TZ = (ri, r2, . . . , r^) and C = (ci, C2, . . . , Cn) are given, where ri > r2 > 
■ . . > and Ci > C2 > . . . > c„. Let the uniquely determined neighbour corresponding to 
{TZ, C) have column sums Vi > V2 > ■ ■ ■ > Vn- Then we define 

1 " 

Note that a(7l,C) is always an integer, since 2a(7l,C) is congruent to 

n n n n 

$^(cj + Vj) = ^Cj + ^Vj = 2^Cj = Q mod 2. 

i=i j=i i=i j=i 

Consider a set F with line sums {TZ,C) and its uniquely determined neighbour Fi. Let 
a = a{TZ, C). It was proved in [31 Lemma 4] that the symmetric difference F A Fi consists 
of a staircases. In this paper we will only use staircases of length 2, which we will define 
below. For the general definition of a staircase, see [3] or 

Definition 2. A staircase of length 2 in F A Fi is a pair of points {pi,P2) in such that 

• pi and P2 are in the same row, 

• pi is an element of F\Fi, 

• p2 is an element of Fi\F. 



3 Main result 

Suppose we are given row sums 71= (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„), 
where ri > r2 > . . . > and Ci > C2 > . . . > c„. Assume that the line sums are consistent 
but do not uniquely determine a set F (hence at least two different sets with these line 
sums exist). Let a = a(7l,C). 

In [3] it was shown that for all F2 and -F3 satisfying these line sums, we have 

IF2 A F3I < Aa^/2\F2\. 

One may wonder how close we can get to achieving this bound. Our theorem shows that 
we can construct two sets that have a symmetric difference of size at least 2a + 2. 



4 



Theorem 1. Let row sums TZ = {ri,r2, . . . ,rm) and column sums C = (ci,C2, . . . ,c„) be 
given, where ri > r2 > . . . > and Ci > C2 > . . . > c„. Assume that the line sums are 
consistent but do not uniquely determine a set F. Let a = a{TZ,C). Then there exist sets 
F2 and F3 with these line sums such that 

IF2 AF3I > 2a + 2. 

This bound is sharp: for each a > 1 there are line sums {TZ,C) with a = a{TZ,C) such that 
for any F2 and F3 satisfying these line sums we have IF2 A F3I <2a + 2. 



4 Proof 

In this entire section, the row sums 71= (ri, r2, . . . , r^) and column sums C = (ci, C2, . . . , c„) 
with ri > r2 > . . . > Tm and Ci > C2 > . . . > c„ are fixed. Furthermore, Fi is the uniquely 
determined neighbour corresponding to (7l,C), and a = a(7l,C). We denote the column 
sums of Fi by t>i > t>2 > . . . > f„. 

The proof is constructive. We will construct F2 and F3 such that they have the desired 
property. We will do this by changing a set F step by step. Only the final result of the 
construction will be called F2 (or F3); the intermediate sets will always be called F or F'. 
In Section |5] the construction is illustrated by an example. 

Let the columns j for which Vj > Cj have indices ji < ^2 ^ • • • ^ ja, where each such j 
occurs Vj — Cj times. Similarly, let the columns i for which Vi < Ci have indices ii < ^2 < 
. ■ ■ < ia, where each such i occurs Cj — Vi times. Define a column pair as a pair {it,jt)- 
The consistency of the given line sums assures that it > jt for all t. For convenience, define 
^0 = Jo = and i^+i = ja+i = n + l. 

We will construct both F2 and F3 by starting from F = Fi and then for each t moving an 
element of F from column jt to column it in the same row. After we have done that for 
t = 1, 2, . . . , a, the row sums of F have not changed, while the columns of F have changed 
from vi, V2, . . . , Vn to Ci, C2, . . . , c„. The symmetric difference Fi A F then consists of 
a staircases of length 2. Each staircase is confined to a single row and corresponds to a 
column pair {it,jt)- We will show that we have a certain freedom in choosing the staircases. 

Suppose we have moved an element for each of the column pairs (^2,^2), 
{it-i,jt-i), where t > 1. The resulting set is called F and has column sums c[, c'2, . . . , c^. 
Now we want to move an element from column jt to column it. For this we need a row / 
such that the point (/, jt) G F and (/, it) ^ F. We have c'j^ > cj^ > ci^ > c[^, so c'j^ > c'^^ + 2. 
Hence there must be at least two rows that contain an element of F in column jt but not 
in column it. This proves the existence of such a row /, and in fact at least two choices for 
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I are possible. Now we move the element {l,jt) to {I, it). The row sums of F do not change, 
while the column sum of column jt decreases by one and the column sum of column if 
increases by one. 

We construct both F2 and F3 using the construction above. First we construct F2, making 
arbitrary choices for the rows in which we move elements. Then we will construct F3. For 
this we let the choices in the construction depend on F2, in a way we will describe below. 

Let Pi, P2, . . . , Pr he the distinct column pairs, where P^ has multiplicity kh'- the column 
pair Pi is equal to each of the pairs {ii,ji), (ifci, jfci), the column pair P2 is equal to each 

of the pairs {iki+i,jki+i), (^fei+fca^ifei+fea); and so on. We have ki + k2-\ h/cr = a. For 

two consecutive column pairs {it,jt) and {it+i,jt+i) that are not equal we have it+i > it, 
jt+i ^ jt or if-^i > if, jt+i > jt, so the second pair contains a column that did not occur in 
any of the previous pairs. This means that in Pi, . . . , at least r + 1 different columns 
are involved. For each Ph, we denote one of the columns in Ph as the final column of Ph in 
the following way. 

• If one of the columns in Ph also occurs in Ph+i, then the other does not occur in 
Ph+i, ■ ■ ■ , Pr- We call the latter the final column of the pair. 

• If both columns in Ph do not occur in Ph+i, ■ ■ ■ , Pr, and one of the columns occurs 
in P^_i, then the other does not occur in Pi, . . . , Ph-i- We call the former the final 
column of the pair. 

• If both columns in Ph do not occur in Pi, . . . , Ph-i nor in Ph+i, ■ ■ ■ , Pr, then we 
arbitrarily pick one of the columns in Ph and call it the final column of the pair. 

By definition, we have the following properties: the final column of Ph does not occur in 
Ph+i, ■ ■ ■ , Pr, and if the other column does not occur in Ph+i, . . . , Pr either, then the latter 
column only occurs in P/j. 

Our goal is to construct P3 in such a way that, for all h, in the final column of P^ the 
symmetric difference between P2 and P3 is at least 2kh, while in any other column that 
occurs in one of the column pairs the symmetric difference between P2 and P3 is at least 
2. (There is at least one such a column, since there are exactly r final columns, while at 
least r + 1 columns are involved in the column pairs.) If we can achieve that, then we have 

IP2 A P3I >2ki + 2k2 + ... + 2kr + 2^2a + 2. 

To achieve this, we choose the rows in which elements are moved for all equal column pairs 
at once. First we choose the rows for all pairs equal to Pi, then for all pairs equal to P2, 
and so on. 
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Let t be the index of the last column pair in a sequence of k equal column pairs 

{h-k+l, jt-k+l) = {it-k+2, jt-k+2) — . . . — {it,jt), 

where {it-k,jt-k) 7^ {it~k+i, jt~k+i) and {it-Jt) 7^ {it+i,jt+i)- Suppose we have moved ele- 
ments already for the cohimn pairs . . . , {it-k,jt-k)- Call the resulting set F, with 
column sums c[, . . . , c'^. Assume that it is the final column of {it,jt) (the case where jt is 
the final column, is analogous). So we have if 7^ if+i- Also, we have one of the following 
two properties: 

(A) Jt = jt+i, 

(B) jt 7^ Jt+i, and jt-k 7^ jt-k+i- 

As this is the last time column it occurs, we need to choose the rows in such a way that 
by moving the elements of F the symmetric difference between F and F2 in this column 
becomes at least 2k. Also, in case (B) we want the symmetric difference in column jt to be 
at least 2. 

Since we need to move k elements out of column jt into column it, we have c'j^ > Cj^ + k > 
C'it + k > c'-^ + 2k, so there are at least 2k rows / such that {l,jt) G F and {I, it) ^ F. Let 
R be the set of those 2k rows. (If there are more than 2k possible rows, then pick 2k of 
them.) We distinguish between two cases. 

Case 1. Suppose there are k different rows / in R such that {I, it) ^ i^2- Then we move 
elements from column jt to column it in each of those k rows. Call the resulting set F' . We 
have {I, it) G F'\F2 for k different values of /. The number of elements of F' in column it 
must be equal to the number of elements of F2 in column it, so there are also k different 
values of I for which {I, it) G F2\F'. Hence the symmetric difference between F' and F2 in 
this column is at least 2k. 

In case (A) we are now done, as column jt will be handled in a later column pair. Suppose 
we are in case (B). The column jt only occurs in the column pairs {it-k+i, jt-k+i), ■ ■ ■ , 
{it,jt), which are all equal. If for a row I we have {I, it) ^ F2, then in the construction of F2 
this row was not used for a staircase corresponding to the column pair {it,jt) (or one of the 
equal ones), so we must have {l,jt) G F2. Hence after moving elements we have k different 
values of / for which {l,jt) G F2\F'. So in column jt the symmetric difference between F' 
and F2 is at least 2k > 2. 

Case 2. Suppose there are at least k + 1 different rows / in R such that (/, it) G F2. Let R' 
be a set of /c + 1 of those rows. Pick one of the rows in R' and call it Iq. Let R" consist of 
Iq and the k — 1 other rows in R\R' (for which it may or may not hold that {I, it) G F2). 
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Move elements from column jt to column it in each of the k rows in R" . Call the resulting 
set F'. Then for all k rows / in R\R" we have {I, it) € F2\F'. Similarly to above, we find 
that the symmetric difference between F' and F2 in column it is at least 2k. 

Again, in case (A) we are done. Suppose we are in case (B). As column jt only occurs 
in the column pairs {it-k+i,jt-k+i), {h,jt), which are all equal, for at most k rows 
Z in i? we have {l,jt) ^ F2. This means that we can choose Iq above in such a way that 
(^O)Jt) £ F2. After moving the elements, we then have {lo,jt) ^ F2\F'. So the symmetric 
difference between F' and F2 in column jt is at least 2. 

At least one of Case 1 and Case 2 above must hold, since there are 2k rows in R. Therefore 
we have finished the construction of F2 and -F3 such that F2 A F3 > 2a + 2. 

We will now prove the second part of Theorem [T] by giving a family of examples for which 
the bound of 2a + 2 is sharp. Let s > 1 be an integer. Take m = n = s + 1 and let all row 
and column sums be equal to 1. These line sums are consistent. The uniquely determined 
neighbour Fi has column sums Vi = s + 1, V2 = = . . . = Vs+i = 0, so a = s. 

Suppose F2 and F3 satisfy the given row and column sums. We have IF2I = IF3I = s + 1, 
hence 

IF2 A F3I < IF2I + IF3I = 2(s + 1) = 2a + 2. 

This completes the proof of Theorem [T] □ 

Remark 1. There do not seem to be very many examples for which the bound of 2a + 2 is 
sharp. In particular, they all seem to have m = n = a + 1. However, even in more general 
cases, when a is much larger than n, the bound is not very far off. Take for example m = n 
and let all line sums be equal to k, where k < ^n. The uniquely determined neighbour has 
k column sums equal to n and n — k column sums equal to Q, so a = k{n — k). As 
n — k > |n, we have a > ^kn. Suppose F2 and F3 satisfy the given row and column sums, 
then IF2I = IF3I = kn, hence 

\F2 A F3I < IF2I + IF3I = 2kn < 4a. 



5 Example 



We illustrate the construction in the proof by an example. Let row sums (5, 5, 5, 4, 4, 2, 1, 1) 
and column sums (6, 6, 6, 3, 3, 3) be given. The uniquely determined neighbour Fi has the 



same row sums, but column sums (8,6,5,5,3,0) (see Figure 1(a)). From this we derive 



that a = 4 and that the four column pairs are (1, 3), (1, 6), (4, 6) and (4, 6). 
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To construct F2, we move one element from col- 
umn 1 to column 3, one element from column 1 
to column 6, and two elements from column 4 
to column 6. We choose the rows to move ele- 
ments in arbitrarily from the available rows. If 
we choose rows 7, 1, 2 and 3 respectively, we ar- 



rive at the set F2 shown in Figure 1(b) 



00000 

00000 

00000 

0000 

0000 

00 

O 

O 



Now we construct the set F3 step- by-step, fol- 
lowing the proof of the theorem. We start with 
Fi, shown again in Figure 2(a) For the first col- 
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(a) The set Fi with 
its row and cohimn 
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6 6 6 3 3 3 
(b) The set F2 with 
its row and column 
sums. 



Figure 1: 



we need to move an element from col- 
umn 1 to column 3. The available rows are 6, 7 
and 8. We need only two of them, so let us take 
R = {7, 8}. Column 3 is the final column in this column pair, so in this column we need to 
make sure that we achieve a symmetric difference of at least 2 with F2. We have (8, 3) ^ F2, 
so we are in case 1 and we pick row 8 for our staircase. Hence we delete the element (8, 1) 
and add the element (8, 3). The new situation is shown in Figure 2(b)[ 
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(a) The set Fi. (b) After one step. 
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(c) After two steps. (d) The set F3. 



Figure 2: The construction of the set F^. 



The next column pair is (1, 6). Now column 1 is the final column pair, and all rows except 



row 8 are available. We are again in case 1 and pick row 4. Figure 2(c) shows the new 
situation, after deleting (4, 1) and adding (4,6). 

Finally, we need to move two elements at once for the column pair (4,6), which occurs 
twice. Column 6 is the final column, so we need to achieve a symmetric difference of at 
least 4 with F2 in this column. We also need a symmetric difference of at least 2 in column 
4 (case (B)). We have R = {1, 2, 3, 5}. As (1, 8), (2, 8) and (3, 8) are all elements of F2, we 
are in case 2. We have R' = {1, 2, 3} and we need to find an Iq G R' such that (/q, 4) G F2. 
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The only possible choice is Iq = 1. We find R" = {1,5}, so we delete (1,4) and (5,4), and 
we add (1,6) and (5,6). This completes the construction of F^. The resulting set is shown 
in Figure [2(d) 



The construction guarantees that the symmetric difference between F2 and F3 is at least 
2a + 2 = 10, but we have in fact constructed two sets with symmetric difference 14. 
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